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Abstract 

The BaBar experiment has completed its first year of data-taking during which 21 fb^ 1 of data 
were accumulated. We present preliminary results on two types of charmless B decays based on 
this data sample. The first, B ^ K*°"f, is a so-called electromagnetic penguin. And the second, 
B° — > h + h ~ , where M •* can be either a pion or kaon, can come from either penguin diagrams or 
Cabibbo suppressed tree diagrams. 
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1 Theory and Motivation 



In the Standard Model, the charmless B decays that will be discussed in this paper proceed via 
the Feynman diagrams shown in Figure |l[ The diagram shown in Figure |l](a) is the "penguin" 
diagram. This diagram has a loop with a u-type quark, which can emit either photon 
(7), or a gluon (g) before recombining to form either d or s quark. Due to its heavy mass, the 
diagram containing a top quark is much larger than those containing u or c quarks, which can thus 
be ignored. Therefore, there are CKM factors of Vt s or Vtd at the recombination vertex, depending 
on whether the final state contains an s or d quark. If the boson emitted in the quark-loop is a 
photon, then final states such as B° — > K* 0n f [] (with an s-quark in the final state) and B° — > p°7 
(with a d-quark in the final state) can be produced. If the emitted boson is a gluon, then fully 
hadronic states such as B°— > K^ir" can be produced. Neither diagrams containing Z° bosons, nor 
those with cf-quarks have yet been observed. 

Measuring decays produced by the penguin diagram is interesting because it can potentially give 
a low energy window on high energy phenomena. For example, in the Standard Model, measuring 
the branching ratio of B° — > p 0r y will give information about the value of Vtd- Beyond the standard 
model, a charged Higgs (H + ) or similar particle could replace the W and either enhance the 
branching ratio for B° — > K*°j or lead to direct CP violation |ffl]. 

The "tree diagram" of Figure |l|(b) can lead to the decay B° — > ir + ir~. This is especially in- 
teresting, because the final state is a CP-eigenstate and the decay has a Matrix Element that is 
proportional to the the CKM factor V u \,. This means that the decay could be used to measure 
time-dependent CP violation. Such a measurement, when combined with constraints from isospin 
conservation, could yield a very clean measurement of the unitarity angle a, which is difficult to 
measure in any other way ||]. However, the penguin diagram can also produce the same final state 
and, since it has a different matrix element, could complicate the measurement. Regardless, a 
precise measurement of a will require higher statistics than are currently available. 

In this paper, we will present results on analyses of the decay modes B° — > isT*°7 and decays 
to two light charged hadrons (B° — > 7r + 7r~, K + ir~, K + K~), which we will often call B°— > h + h ~. 
Preliminary BABAR results on charmless decays based on about half of the current data set have 
been presented in (B° — > h + h _ ) and in Q (B° — > K* ^). Here, we present updated preliminary 
results based on the full data-set of the first year. 

2 The PEP-II B Factory 

The data that will be described in this paper was taken with the BABAR detector at the PEP-II B 
Factory, which is at located the Stanford Linear Accelerator Center (SLAC) and has been running 
since May of 1999. PEP-II is an asymmetric e + e~ collider that collides a beam of e~'s with an 
energy 9.0 GeV with a beam of e + 's with an energy of 3.1 GeV. This produces a center of mass 
energy of equal to the T(4s) mass and thus provides a copious source of B mesons. The peak 
luminosity achieved in "Run 1" , which stretched from November of 1999 to October of 2000, was 
3.3 x 10 33 cm _2 sec— 1. The integrated luminosity collected during that time was approximately 21 
/ft -1 , which is equivalent to (22.7 ±0.4) x 10 6 BB pairs. Additionally, 3/6 1 of data was taken at 
energies below the T(4s) peak. 

The asymmetric energies of the PEP-II beams are required for doing the time-dependent 
CP asymmetries that are the heart of the BABAR physics program ||. However, they cause 

1 Charge conjugation is assumed throughout this paper unless specifically stated. 




Figure 1: Feynman diagrams responsible for charmless B decays, (a) is the "penguin" diagram, 
which contributes to the all of the final states being discussed in this paper, (b) is the "tree" diagram, 
which can contribute to the B° — » 7t + 7t~ decay. 

some problems for the measurement of charmless branching ratios. For example, the range of 
momenta of the daughters produced in two-body B decays is rather narrow in the center of 
mass frame (2.4 GeV < P2-body,CMS < 2.8 GeV). In the lab frame, however, the range is wider 
(2GeV < P2body,iab < 4 GeV). This means, for example, that the the particle identification device 
used to separate 7r's from if's needs to work over a wider range of momenta than it would otherwise. 

3 The BABAR Detector 

Figure ^ shows a schematic of the BABAR detector. The innermost component of BABAR is 
5 layer Silicon Vertex Tracker. Outside of that is a 40 layer central Drift Chamber. The DIRC 
(Detector of Internally Reflected Cherenkov light) ||, which is used for charged hadron identifi- 
cation is composed of two parts: the quartz radiator bars, which are located just outside of the 
Drift Chamber, and the array of Photo Multiplier tubes, which is located on the backward side 
of detector outside the central region. Surrounding the DIRC radiator bars is an electromagnetic 
calorimeter using Csl crystals. Outside of this is the Superconducting Solenoid, which provides a 
1.5 T field. Finally, the solenoid is surrounded by an iron flux return, which has gaps instrumented 
with resistive plate chambers to provide muon and neutral hadron identification. The capabilities 
of BABAR are described in more detail in ||. Since the charged hadron identification provided by 
the DIRC is key to both of the analyses described in this note, we will describe it in more detail 
here. 

4 The DIRC 

Figure |3] shows the principle of DIRC operation. As charged particles pass through the radiator 
bars, Cherenkov light is produced. A fraction of that light is trapped inside the rectangular bars 
by total internal reflection. As the light than propagates down the length of the bar, typically 
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Figure 2: Side view of the BABAR detector. 
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Figure 3: The principle of DIRC operation. See text for explanation. 

bouncing 200 times, the Cherenkov pattern is preserved, up to reflection ambiguities. At the end of 
the bar, the light emerges into a stand-off region where the Cherenkov pattern is allowed to expand 
before being detected by an array of roughly 11,000 phototubes. 

Measurement of the Cherenkov angle of the detected light provides excellent w, K separation. 
The performance of the DIRC can be checked using samples of 7r's and K J s identified in D* + — > 
ir + D°(D — > K~ir + ). Figure §(a) shows the performance of a K selection algorithm based on the 
DIRC Cherenkov angle measurement on these samples. The efficiency for correctly identifying ICs 
is roughly 90%, while the probability of mis-identifying a n as K is roughly 2 % at low momentum, 
but degrades somewhat at high (> 2.5GeV) momentum. Figure |4](b) shows the bands of Cherenkov 
angles measured by the DIRC for high momentum tracks. The bands grow closer together at high 
momentum, but the separation is still better than 2.5 a at 4 GeV, which is the kinematic limit for 
two-body B decays. 



5 Analysis Procedure 

A very similar procedure is used for both the B° —* K*°j and B° — > h + h'~ analyses. The first 
step is to compose "B candidates" of the desired final state out of the available charged tracks, 
photons and 7r°'s. Then, background is rejected based on event shape variables. Background levels 
are measured using off-resonance data or sideband regions of on-resonance data. To determine the 
number of signal events, a fit is performed to a set of kinematic variables. The signal efficiency is 
calculated in Monte Carlo, including adjustments for tracking efficiency and other related effects. 




Figure 4: ir/K separation of the DIRC. a) shows the performance of a kaon selection algorithm 
using DIRC information, b) shows the raw Cherenkov angle measurement for samples of it's and 
K's. 

6 Background Suppression 

Because the cross-section for continuum light-quark production is much higher than that for BB 
production cr(e + e~ — ► qq) « 3o"(e + e~ — > BB) and because the tracks coming from continuum events 
tend to have higher momentum, the background to charmless decays comes mostly from continuum 
events. This background may be suppressed by noting that in the center of mass frame, continuum 
events are more "jetty" than the more spherical BB events. 

The first step in this procedure is to define the "axis" of the candidate B-decay. In the case 
of B° — » if*°7, this is taken to be the direction of the photon. While for the case of B° — ► h + h ~, 
it is taken as the thrust axis of the B candidate. Then, the tracks and neutral energy of the the 
remainder of the event (excluding the B candidate) are examined. In the case of B° — > K* ^/ the 
thrust axis of the remainder is calculated and the angle between it and the candidate axis is called 
Of. Figure ||(a) shows the distribution of the cos Ot for signal and background. As expected for 
the "jetty" background, cosOt is peaked toward 1. For the more spherical BB events, it has a flat 
distribution. 

In the case of B° — > h + h~ a more sophisticated approach, which was originally developed by 
CLEO is used 0. This method defines a set of 9 concentric cones centered on the candidate axis. 
Each cone subtends 10° of solid angle and is folded to combine the forward and backward intervals. 
The energy in the ith cone is called x% and a "Fisher Discriminant" , F, is defined as 

9 

F = ^2aiXi. (1) 

i=l 

The coefficients a, are varied so as to maximize the separation between signal and background. 
Figure ||(b) shows histograms of this quantity for signal and background. 
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Figure 5: Event shape variables used for background suppression, a) shows the cosOt variable used 
in B° —> K*°~f. b) shows the Fisher Discriminant used in B°—> h + h~ . Solid circles are for signal, 
open circles for background. 

7 Kinematic Variables 

Two kinematic variables are traditionally used in exclusive ^-reconstruction in symmetric colliders 
operating at the T(4s). The first, known as the "beam-constrained mass" is defined as, 

m B = jELm-PB, ( 2 ) 

where E oeam is the beam energy and ps is the momentum of the B candidate in the center of mass 
system. Since the beam energy is precisely known, this variable has much better mass resolution 
than would be achieved by using the measured energy of the B candidate. The second variable is 
defined as 

AE = E B — E beam , (3) 

where Eb is the measured energy of the B candidate. In the case of B ^ h+h'-, if one assumes 
that both daughters have the pion mass, then AE will peak at zero for signal events for which both 
daughters are pions. If one or more of the daughters are not pions, then the peak will be shifted 
away from zero, thus providing particle identification. In the case of B° — ► K*°j, where there is no 
ambiguity in the particle masses, one expects a peak at zero. 

In an asymmetric collider, one can simply boost the track parameters back into the center of 
mass frame and calculate as usual. However, this requires assignment of masses to each of 
the particles. In some analyses, such as B° — ► h + h ~, this is undesirable and can lead to worse 
m,B resolution if the wrong mass assignment is made. In such cases, it is desirable to define a new 
variable called the energy-substituted mass mgs, 



mES = \j(-s + p -pB)/E$-p B , (4) 



where y/s is the center of mass energy, p~o and Eq are the three-momentum and energy of the T(As) 
and p~B is the three-momentum of the B candidate. This variable has the advantage that it uses 



Figure 6: K*° mass peak for the final B° — > K* 0/ j sample. 

only quantities measured in the lab system, and thus requires no mass assignments to be made. It 
is identical to if the lab and center of mass frames coincide. It is also equivalent to if there 
is no ambiguity in the particle masses. A new version of AE is also defined, 

AE* = E* B - (5) 
which also has shifts for different daughter masses. 

8 B ^ K*°i Analysis 
8.1 B ^ K*°~f Candidates 

The B° — > K*°j decay is searched for in the mode K*° —> K n^. B° — > K*°j decay candidates are 
composed out of available photons and K* 0, s, which are in turn composed of K s and tv^'s. 

In order to be considered, the photon candidate must be located cleanly within the active region 
of the calorimeter (—0.73 < cos6i a b < 0.9) and must have an energy close to that expected in the 
center of mass (2.3 GeV < E lt cMS < 2.8 GeV). In addition, the photon must not be consistent 
with having come from a 7r° or r\ decay. Two cuts are applied to insure this. The transverse shape 
of the shower is used to reject those 7r°'s for which the two photons are not spatially separated. 
And, to reject 7r°'s and ry's that do form two separate calorimeter clusters, the candidate photon is 
matched up with all other photons in the event and rejected if any of them gives a mass consistent 
with m n o or m„. The efficiency of these photons cuts is 77%. 

Similarly, the K*° candidate must pass a number of cuts. First, the K + and tt~ candidates must 
pass DIRC particle identification cuts. Then, the mass of the Kir system must be consistent with 
the K*° (0.796GeV < vtlk-k < 0.996GeV). Finally, the helicity angle, defined as the angle between 
the K + direction in the K* rest frame and the K* flight direction, must satisfy | cos 0/^1 < 0.75. 
The efficiency of the K*° cuts is 56 %. Figure |6| shows the K*° mass peak in the final B° — > K* 0r y 
sample. 
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Figure 7: The final sample of the B° — > K*°j search, a) is a scatter plot of AE vs. tub- b) is a 
projection onto the tub axis, showing a clear B° — > K* 0r y signal. 

8.2 B -> ir*° 7 Mass Fit 

Once the photon and if* candidates have been identified, B candidates are formed from them. 
These candidates must pass a number of cuts in order to be considered in the final sample. The 
AE* of the candidate, defined in Equation |] must pass a cut of -200MeV < AE* < lOOMeV. 
The "thrust angle" must pass the cut | cos 6t\ < 0.8. And finally, the polar angle of the B candidate 
in the center of mass frame must pass the cut | cos#|j| < 0.75. The events that pass these cuts are 
shown in Figure 0(a) and (b). The concentration of events near AE* = and Mb = 5.28 in Figure 
0(a) indicates a strong B° — » K*°j signal. In order to determine the number of signal events, a fit 
is performed on the Mb distribution. As shown in Figure @(b), he background is modeled with an 
"Argus" function, whose shape is determined using data taken off the T(4s) resonance. The shape 
of the signal is Gaussian, and the mean, a and amplitude are floated in the fit. The number of 
signal B° — > K* 0r y events extracted from the fit is N s i g = 139.2 ± 13.1. 

8.3 B° -> K*°~f Results 

The efficiency for observing a B° — ► K* ^/ event is calculated in Monte Carlo. This includes correc- 
tions for known differences between data and Monte Carlo for tracking efficiency, photon detection 
efficiency and particle identification efficiency. The efficiency is found to be e = 0.209 ± 0.013, 
where the error is purely systematic. One can then use this efficiency to calculate the branching 
ratio as: 

BR(B^K*^)= "*> , (6) 

^BB x 6 x B K* 

where N B g is the number of BB pairs in the sample, and Bk* is the branching ratio for K*° — > 
K~tt + . The systematic error comes mostly from data-derived efficiency corrections. The prelimi- 
nary result for the branching ratio is then: 

BR{B° - K*° 7 ) = (4.39 ± 0Al stat ± 0.27 syst ) x lO" 5 (7) 



Measurement of a possible CP asymmetry in this decays is also of interest. The asymmetry is 
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Figure 8: Signals for three other modes of B^ K*"f. 



defined as: 



N(B° —* K* Qr y) — N(B° — >■ if *°7) 



(8) 



Many of the systematic errors of the branching ratio measurement cancel out in this measurement 
and we are left with a preliminary measurement of: 



A C p 



-0.035 ± 0.094 ± 0.022 



(9) 



8.4 Other B° -> K* 7 Modes 



Figure || shows evidence for three other modes of B — ► K*y. B° — > if*°7,if*° —y 7r°K®; B + — > 
if* + 7, K* + — > if + 7r°; £? + — ► if* + 7, if* + — > if|?7r + . Although clear signals are evident in each of 
these modes, their efficiencies and systematic errors have not yet been fully evaluated. For this 
reason, no branching ratios or asymmetries will be presented for them at this time. 



Variable 



Definition 



Energy substituted B Mass, described in Section [7] 
Energy difference, described in Section |7| 

Fisher discriminant output, described in Section ^ 
DIRC Cherenkov Angle for positive track 
DIRC Cherenkov Angle for negative track 



TIES 

AE 
F 

®c,+ 



Table 1: Variables used in the Maximum Likelihood fit of B° — > h + h' and their definitions. 

9 Analysis of B ^ h + h'~ 

9.1 B ^ h+h! Candidates 

B° — > h+h'~ candidates are composed out of pairs of oppositely signed charged tracks, both of 
which must pass tracking quality cuts and DIRC quality cuts. The candidates must pass a cut 
on the "sphericity angle" which is almost identical to the "thrust angle" , 6t defined in Section |(| 
of cos 6 S < 0.9. Furthermore, the candidate events must pass cuts on the second Fox- Wolfram 
moment of R2 < 0.95 and on sphericity of S > 0.01. Finally, loose cuts are placed on the kinematic 
variables of the decay: 5.2 GeV < m ES < 5.3 GeV and -0.15 GeV < AE < 0.15 GeV. In Monte 
Carlo, these cuts are found to have an efficiency of e ~ 0.45 for the B°^> 7r + 7r _ mode. The 26404 
events that pass these cuts in the data sample are put into a Maximum Likelihood fit described in 
the next section. 

9.2 B ^ h + h' Likelihood Fit 

Since a small number of signal events are expected in the B°— ► h+h "modes, it is desirable to use a 
Maximum Likelihood fitting technique, rather than the more traditional "cut-and-count" technique. 
For the purposes of the fit, each event is considered to come from one of eight hypotheses: a true 
B° — > ir+ir~ event, a true B° — > K+tt~ event (or charge conjugate), a true £>°— > K + K~ event, 
a background B°^tt + tt~ event, a background B°^K + it~ event (or charge conjugate), or a 
background B° — > K + K~ event. The parameters that are varied in the the likelihood fit are the 
number of events in each of these hypotheses. The event variables that are fit to are listed in Table 



The likelihood function is defined as follows. A probability distribution function (PDF) is 
calculated for each event and hypothesis. It has the form: 

phypo _ phypo phypo phypo phypo phypo 
r event ~ r M ES 'AE F r e c ,+ 9c,- ' 

where the P's are the PDF's for each variable given an event hypothesis. Typically, these are 
parameterized in terms of simple functions. The full likelihood function is then written, 

N 

C = exp(- J2 Nhypo) II [ E N hypo P^Z] , (11) 

hypo j=l hypo 

where Nh ypo is the number of events in each hypothesis and N is the total nmber of events in the 
sample. This function is then maximized using standard tools. 



(10) 
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Figure 9: tyies distributions used to calibrate the ities PDF's. a) is from B~ — > D (— > i"C~7r + )7r~ 
cfaia and is used /or i/ie signal PDF. b) is from AE sideband data and is used for the background 
PDF. 

9.3 Calibration of PDF's 

In order for the the fit to be valid, it is essential that the individual PDF's properly describe 
the distributions for the signal and background hypotheses. These PDF's are calculated based on 
calibration samples of real data and checked with Monte Carlo simulation. The following sections 
describe this procedure in more detail for each of the five fit variables. 

9.3.1 m E s PDF's 

The shape of the uies distribution for signal is taken to be Gaussian. Since its width is dominated 
by the beam spread energy, it can be reliably calculated based on fully reconstructed decays of the 
type B~ —* > K~ir + )7r~ . Figure ||(a) shows the shape of this distribution. For background, 
the shape is assumed to be the "Argus Function", which is fitted to data taken on the T(4s) peak, 
but with AE outside of the signal region. Figure ^(b) shows this distribution and its fit. 

9.3.2 AE PDF's 

The AE shape for signal is assumed to be Gaussian. In contrast to the tues case, however, the 
width is dominated by tracking resolution and is therefore different for B~ — > Z)°(— > K^tt + )-k^ than 
for B ^ h + h'~. In this case, we compare the width measured in data for B — > D°(— > K 7r + )7r , 
o~ae = 19MeV, with that found in Monte Carlo, a^E = 15MeV. We then scale up the width 
found in signal Monte Carlo (o~ae = 21MeV) by the same factor to obtain the width used for the 
PDF, o~t\E = 26 ± 5MeV. For background, the shape is taken to be a polynomial, which is fitted 
to data taken on peak, but outside of the AE signal region. 



Decay Mode 


^signal ± & stat ± 0~syst 


BR BABAR (xlO" 6 ) 


K+TT- 

K+K- 


41 ± 10 ± 7 
169 ± 17±i| 
8.2i^ ± 3.3 


4.1 ± 1.0 ±0.7 
16.7 ± 1.61:? , 7 
< 2.5(90% C.L.) 



Table 2: Number of events for each mode found in the Maximum Likelihood Fit and the corre- 
sponding measured branching ratios. 



Parameter 


Measurement 


BR( J B U ^ iT* u 7) 
Acp{B°^K*°j) 
BR^ ^ K+7T ) 

BR(i? ^ 7T+7T-) 


(4.39 ±0.41 ±0.27) x 10~ 5 
-0.035 ± 0.094 ± 0.022 

(16.7 ± l.etla) x 10~ 6 
(4.1 ±1.0 ±0.7) x 10~ 6 



Table 3: Summary of measurement presented in this paper. In each case, the first error is statistical 
and the second s systematic. 

9.3.3 Fisher Discriminant PDF's 

Double Gaussians are used for the Fisher discriminant shape. Signal Monte Carlo is used to 
determine the parameters for signal. It is checked with the B~ — > D (— ► K~tt + )tt~ sample. For 
background, the parameters are determined using the tues sideband and are checked using off- 
resonance data and continuum Monte Carlo. 

9.3.4 9 C PDF's 

A sample of D* + — ► tt^ 1ow D°(D° — > K ~ir + ) was used to provide clean samples of 7r + 's and K +1 s. 
The DIRC 9c response was parameterized by a Gaussian with mean and width dependent on dip 
angle. Also, it was necessary to include small non-Gaussian "satellite" peaks in order to adequately 
model the performance. 

9.4 B ^ h + h'- Fit Result 

The maximum likelihood fit yields the results shown in Table ||. The systematic errors are calculated 
by varying the parameters of the PDF's, both within their statistical errors and to cover any 
disagreements between data and Monte Carlo. More information about this analysis may be found 
in §. 

10 Conclusions 

Based on a first year sample of 22.4 x 10 6 BB pairs, BABAR has the preliminary measurements 
shown in Table ||. As more data is collected and more decays modes analyzed, BABAR will have 
many more results on charmless B decays. 
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